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(N 1 ABSTRACT 

Gravitational lensing causes a correlation between a population of foreground large- 
' scale structures and the observed number density of the background distant galaxies 

as a consequence of the flux magnification and the lensing area distortion. This cor- 
relation has not been taken into account in calculations of the theoretical predictions 
of the cosmic shear statistics but may cause a systematic error in a cosmic shear 
measurement. We examine its impact on the cosmic shear statistics using the semi- 
, analytic approach. We find that the lensing magnification has no practical influence on 

J> ■ the cosmic shear variance. Exploring possible shapes of redshift distribution of source 

galaxies, we find that the relative amplitude of the effect on the convergence skewness 
is 3% at most. 
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1 INTRODUCTION 



45 

q . The cosmic shear — coherent distortions in distant galaxy images due to the weak lensing by large-scale structures — is now 
' recognized as a powerful tool to measure the mass distribution in the universe as well as a promising way to measure the 
CZ2 . cosmological parameters (Mellier 1999; Bartelmann & Schneider 2001 for reviews). Although its signal is very weak, recent 
reports on the detections demonstrate that a well developed data analysis algorithm has been established (Van Waerbeke et 
al. 2000; Wittman et al. 2000; Bacon, Refregier & Ellis 2000; Kaiser, Wilson & Luppino 2000; Maoli et al. 2001). 

So far, the detections were obtained from relatively small fields, which limit the statistical analyses of the surveys to second 
order statistics (the variance or two-point correlation function of the cosmic shear). Their amplitude reflects the amplitude of 
' the density contrast, and thus provides a constraint on the combination of the values of f2 m and erg (Bernardeau, Van Waerbeke 
& Mellier 1997, hereafter BvWM97; Jain & Seljak 1997, Maoli et al. 2001). The skewness of the lensing convergence is, on 
the other hand, known to be sensitive to f2 m almost independently on erg (BvWM97; Van Waerbeke et al. 2001b, hereafter 
vWHSCBOl), and thus it may break the degeneracy among the cosmological parameters constrained from the second order 
statistics. A precise measurement of the skewness is, therefore, one of main goals of on-going wide field cosmic shear surveys 
such as the DESCART project^. 

It was pointed out by Hu & Tegmark (1999) that such wide field cosmic shear surveys have a potential for probing the 
cosmological models as accurately as the cosmic microwave background. In order to get their ability to the fullest, there 
are, however, some issues which should be developed/examined in details: (i) making an accurate theoretical prediction of 
the cosmic shear statistics at intermediate scales (0.5-10 arcmin). On such scales the signals will be detected easily, however 
neither perturbation theory nor the hierarchical ansatz apply (e.g., Jain & Seljak 1997; vWHSCBOl). (ii) examining possible 
corrections which arise from higher order correction terms in calculations of the theoretical predictions, for example, Born 
approximation and lens-lens couplings (BvWM97; Schneider et al. 1998; vWHSCBOl), and the source clustering (Bernardeau 
1998; Hamana et al. 2000). They are especially important for the convergence skewness because it is (in the perturbation 
theory sense) a quantity of 4th order of 5^ (see §2 for details), (iii) examining the impact of the intrinsic shape correlation of 
source galaxies; in the theoretical calculations of the cosmic shear statistics, it is supposed to be negligible (Croft & Metzler 
2000; Heavens, Refregier & Heymans 2000; Catelan, Kamionkowski & Blandford 2001; Crittenden et al. 2000a; 2000b). (iv) 



* For more information about DESCART project, see http://terapix.iap.fr/Descart/ 
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developing a robust way to correct defects in instruments, in particular the point spread function anisotropy (Kaiser Squires 
& Broadhurst 1995; Kuijken 1999; Erben et al. 2001; Bacon et al. 2000). 

In this paper, we focus on the lensing magnification effects on the cosmic shear statistics, which have not been taken 
into consideration in the theoretical calculation of the cosmic shear statistics but may cause a systematic error in their 
measurements. The lensing magnification has two effects on a deep galaxy survey: One is the flux magnification; the lensing 
changes the apparent galaxy size, leaving the surface brightness invariant, therefore the flux of a distant galaxy is magnified^]. 
The other is the area distortion; the lensing also changes the unit solid angle at the source plane, and thus the number density 
of the distant galaxies varies with the direction in the sky even if the intrinsic source distribution is uniform. As a consequence, 
the lines of sight to distant sheared galaxies may not be random lines of sight in presence of the lensing magnification. In 
fact, the lensing magnification effects are prominent within the galaxy cluster region where the number density of distant 
galaxies in the optical or near IR bands is measured to be smaller than the average value. Furthermore, the variation of the 
galaxy number density as a function of the distance from the cluster center (so-called, depletion curve) were measured in 
some distant clusters of galaxies (e.g., Broadhurst, Taylor & Peacock 1995; Fort, Mellier & Dantel-Fort 1997; see also Mellier 
1999 for a review). These observational facts indicate that the lensing magnification induces an (anti-)correlation between the 
distribution of the distant galaxies in the sky and the population of the lensing structures, i.e., the number density of the 
galaxies behind a lensing structure tends to be smaller than the average number density, whereas that behind a void tends to 
be larger than the average value. 

The purpose of this paper is to quantitatively examine the lensing magnification effects on the cosmic shear statistics, 
especially on the convergence skewness. To do this, we use the nonlinear semi-analytic approach, i.e., the perturbation theory 
approach combined with the nonlinear fitting formula of the density power spectrum (Jain & Seljak 1997; Hamana et al. 2000; 
vWHSCBOl). We focus on correction terms which arise from the presence of the lensing magnification, and are not concerned 
with other correction terms due to, e.g., the lens-lens coupling (vWHSCBOl) and the source clustering (Hamana et al. 2000). 

The outline of this paper is as follows: The calculations of the moments of the lensing convergence in presence of the 
lensing magnification are made in §2. In §3, the effect on the convergence skewness is quantitatively examined in three cold 
dark matter (CDM) models with realistic models of the redshift distributions of the source galaxy. We conclude in §4. 



2 THE SEMI-ANALYTIC APPROACH 

2.1 Fluctuation in a galaxy number count due to the lensing magnification 

Let n s (> S, z) be the unlensed number density of galaxies with redshift within Az of z and with flux larger than S. Then, at 
an angular position </> where the lensing magnification is fi(z, </>), the number counts are changed by the lensing magnification 
effects as (e.g., Bartelmann & Schneider 2001), 

n° bs (> S,z,4>) = f 1 ,, n s (> . S ,, ,z) . (1) 

Supposing that the number counts of the faint galaxies can be approximated by a power-law over a wide range of fluxes i.e., 
n s {> S, z) = n(z)S~ a ( z \ the lensed number counts are rewritten as, 

n° bs (> S, z, 0) = n s (> S, z) n(z, ^J-W- 1 . (2) 

The magnification factor relates to the lensing convergence (k) and shear (7) by = (1 — re) 2 — 7 2 |. We now rewrite 
the lensed source counts as n° bs (> S,z,<j>) = n s (> S,z)[l + Sn s (> S,z, (/>)]. Taking advantage of the power-law form of the 
number counts and also taking the weak lensing approximation (k<1, 7< 1), the fluctuation in the number counts due to 
the lensing magnification is given by, 

8n s {z,(j>) = (i(z, </>) a(z)_1 — 1 

~ 2[a(x)-l] K .(z,4>), (3) 

where n s (z,(j>) is the lensing convergence at an angular position <f> for a source with redshift z (Mellier 1999; Bartelmann & 
Schneider 2001), 



3f2 m Ho 



K s(z, cj)) = — / dxg{x,Xs)$(x,<l>) , (4) 



t We should here note that the lensing causes not only magnification (jj, > 1 with denoting the magnification factor) but also 
demagnification (fi < 1). Throughout this paper, following the usual convention, we use the term "magnification" irrespective of the 
value of the magnification factor. 
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where g is the so-called lensing efficiency function defined by, 

/ X #0 fix) f(X* ~ X) /jr\ 

g(X ' Xs) = ~ fiXs)aix) ■ (5) 
Here x denotes the radial comoving distance, a is the scale factor normalized by its present value, and fix) denotes the 
comoving angular diameter distance, defined as fix) = K~ x ^ 2 sin K 1 ^ 2 x, Xi i~ K)^ 1 ^ 2 sinh(— K) 1 ^ 2 x for K > 0, K — 0, 
K < 0, respectively, where 7f is the curvature which can be expressed as K = (_ffo/c) 2 (fi m + S1a — 1)- The lensing convergence 
is therefore the projected density contrast weighted by the distance combination and the scale factor along the line of sight 
to a source. Note that, in the expression (||), the fluctuation does not depend on the flux because of the power-law form of 
the number counts. In the following sections, we will therefore not explicitly denote the flux dependence of the source counts. 



2.2 Cosmic shear statistics in the presence of the lensing magnification 

Let us consider the measured convergence that results from averages made over many distant galaxies located at different 
distances. Denoting the smoothing scale by 9, such an average can formally be written as, 

Ke = ^iv 3 ... 77 : ' w 

E<=i We Oi) 

where We{x) denotes the weight function of the average, N s is the number of source galaxies, and <f> t and Zj are the direction 
and redshift of i-th source, respectively. For the weight function, the angular top-hat filter (BvWM97) and/or compensated 
filter (Schneider et al. 1998) are commonly adopted (e.g., Van Waerbeke et al. 2001b). In what follows, we consider the top-hat 
filter for the weight function, and in this case equation ([]) is reduced to Kg — Ei=?i K s(zi, (fr^/Ni , where iVf is the number 
of source galaxies within an aperture 6 centered on a direction . . Taking the continuous limit for the source distribution^ , 
equation (jfjj) can be rewritten by, 

J d 2 <p We (0) /*" d X k s ( X , 0K bs ( X , 0) 

fd^Wei*) $™ d X n?°ix,4>) ' () 

where xh is the distance to the horizon and n° bs (Xi 4>) ls > the source number count defined by equation ^ which effectively 
describes the redshift distribution of the source. 

In what follows, the distance distribution of the unlensed number counts is supposed to be normalized to unity, J H dx n s (x) 
1. Let us now expand equation (jjj) in terms of S using the perturbation theory approach (BvWM97). The presence of the 
lensing magnification does not change the expression of the first order term, 

kP = [ d 2 ^W e {ct>) r d Xn 3 ix) r d X9ix',x)5 {1 Hx,<t>)- (8) 



2c 

The second order convergence consists of two terms: One comes from the second order density perturbation, which is formally 
written by replacing the subscript ^ in the first order expression (^) with ^ (BvWM97). The other is due to the lensing 
magnification, 



mag. (2) 



2 £3tWfo^ a I d 2^ Wg((j)) r ff dxns(x)[Q(x) _i] rdx'g(x\x)S {1) ix',<P) f X d x "g(x",x)S (1) (x"A) 

-4 1 — / d 2 4>Wei<t>) d X n s ( X )Hx) - 1} ^ d X ' g( X ' , x)S il) (x'^) ■ (9) 



Using the small angle approximation (Kaiser 1992), equation (|8|) is rewritten in terms of the Fourier transform of the density 
contrast, S(k), by, 

_ 30^ f XH dxns{x) [ X d x 'gix',x) [ ^Wlfix')k ± e]5^[k-, X ']exp{ik x fi X ')], (10) 



2c J — J Q - (2. 

where the wave vector k is decomposed into the line-of-sight component k x and perpendicular to it, fej_, and W(x) is the 
Fourier transform of the weight function. In the case of the top-hat filter, W(x) = 2J\(x)/x where Ji is the spherical Bessel 
function. In the same manner, equation (^) is rewritten by, 

^.(2) = 2 ^3f^Ho^ r H dxn s (x)Kx)-l] f d^gix ,X) f d X " gix" ,X) 



•t See Bernardeau (1998) for a discussion on this point. 
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^5 (1 >[fc'; X ']exp[a;' x /(x')] / [k" ■ X "] e X p[ik" x f( x ")]W[\n X ')k'± + f(x")k"± \6] 



{2nf l " >a.J"^l"»x.<vwj j (2yr 
-^^f 1 £" d X n s (x)[a{x)-l} f d X ' g(x',x) J -0^W{f( X ')k' ± 8]S^[k'; X '] exp[ifc' x /( x ')] • (H) 

The lensing magnification effect makes the average of the convergence non-zero, 
I* \ — /„ ma s-( 2 )\ 



jf d X nsM[a( X )-l]J o d x 'g\ X ', X )M X ',0) 
2(^^) 2 r d X n s ( X ) d X 'n s ( X ')[a( X ')-l] F d X " g( X ", x )g( x " , X ')I ( X ",8) , (12) 



where 

Ux,0) = ^ / dftfcW 2 [/(x")^]^in(x,fc) , (13) 

with Pn n (x,fc) being the linear matter power spectrum. The variance is not affected by the lensing magnification and is given 
by, 

V„{fi) = - M) 2 ) = (4 1)2 ) = (^^y J*" d X n s ( X )J\ X 'g 2 ( X ', X )I„( X ',9). (14) 

In presence of the lensing magnification, the skewness parameter, defined by S3 (6) = {(Kg — (Kg)) 3 ) /V 2 (0), consists of two 
terms: One comes from the second order perturbation, (Kg) 2FT = 3{Kg Kg ) (see BvWM97 and Hamana et al. 2000 for the 
explicit expression). The other arises from the lensing magnification, 

i4r as = <(4 15 + *r M2) - (^)) 3 ) = 3<4 X) V s ' (2) > - 3(4 i)2 >( K r s - (2) > 

i^^y r a d X n siX )H X ) - 1] 

•/ 

-1/ K («)xl2(^) 2 d X n s { x ) ^ dx'n.(x')[«(x')-l] d X " g( x " , X )g( X " , X ')I ( X " ,8) . (15) 



2 
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d X 'n s ( X ) / dx"ff(x",x)5'(x",x') J o(x",0) 



To derive the last expression, we used an approximation, (Msin^Vydfe + fc'|) ~ 27rW / (fc)W / (fc') with 1? being the angle 
between the wave vectors k and k' . For the top- hat window function, the error this approximation induces is extremely 
weak, for instance it is less than 1% for n ~ —1.5 with n being the effective power-law index of the matter power spectrum 
(Bernardeau 1998). Notice that, in the case of a( X ) — 0, which corresponds to the case that sources are selected on a surface 
brightness criterion (see §4 for discussion on this point), the second term reduces to 12V 2 (9). As this immediately suggests, 
the cosmology dependence in the skewness correction term is weak, actually f2 m in the coefficients of equation ( |l5| ) is canceled 
out by that in V 2 (6) (this point is demonstrated in Figure [f]). 

The above calculations are based on the perturbation theory approach. It is well known that on sub-degree scales the 
nonlinearity in the evolution of the density field is very important for the cosmic shear statistics (Jain & Seljak 1997; 
wVHSCB). We take the nonlinearity into account in our computations in the following way: For the variance, the effect of 
nonlinear evolution of the density power spectrum can be included by replacing the linear power spectrum with the nonlinear 
power spectrum, i.e., Py ln (a, k) — > Pnl(i, k) (Jain & Seljak 1997). We use the fitting formula for the nonlinear power spectrum 
given by Peacock and Dodds (1996). For the skewness correction term, all density contrast terms needed for its calculation, 
equation (|ll]) , correspond to the linear order. This is the same situation as for the variance. Following the procedure used for 
this latter case, we simply replace the linear power spectrum with the nonlinear one to include nonlinear effects. We adopt 
the semi-analytic calculation of the skewness in the nonlinear regime developed by vWHSCBOl, which is based on the fitting 
formula of the density bispectrum by Scoccimarro & Couchman (2000). 

It should be here noted that, comparing equation ( p^ ) with (|l4|), the amplitude of the shift in average of the convergence 
from zero is found to be the same order of the variance, and thus is of order of O(10 -4 ). This shift has no practical effect on 
cosmic shear statistics because the constant shift has no effect on the second and higher order statistics by definition. 



3 NUMERICAL RESULTS 

In this section, we numerically examine the lensing magnification effect on the convergence skewness. We take three CDM 
models, a flat model with (ACDM) and without cosmological constant (SCDM) and an open model (OCDM). We use the 
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Table 1. Cosmological parameters. 



Model 


S7 m 


n x 


h 


CT8 


SCDM 


1.0 


0.0 


0.5 


0.6 


OCDM 


0.3 


0.0 


0.7 


0.85 


ACDM 


0.3 


0.7 


0.7 


0.9 



Table 2. Parameters in models of the redshift distribution of source galaxies (a, b, z») and their characteristics (the mean redshift (z) 
and the root-mean-square of the distribution Az) 



Model 


a 


b 


z* 


<*> 


Az 


A 


2 


1.5 


0.798 


1.2 


0.572 


B 


3 


1.8 


0.813 


1.2 


0.456 


C 


5 


2.5 


1.11 


1.5 


0.400 


D 


2 


1.5 


0.598 


0.9 


0.429 



fitting formula of the CDM power spectrum given by Bond & Efstathiou (1984) normalized by the local galaxy cluster 
abundance (Eke, Cole & Frenk 1996; Kitayama & Suto 1997). The parameters in the models are listed in Table 
We assume that n s (z) takes the form, 

Mz) = ,.r[(i + a)/fr] (£r exp 

where T(x) is the Gamma function. We explore four models for the shape of the distribution. The parameters as well as 
characteristics (the mean redshift (z) and the root-mean-square of the distribution Az) in each model are listed in Table [| 
Note that only model A matches roughly the observed redshift distribution of galaxies in current cosmic shear detections 
(Van Waerbeke et al. 2000). However, we test the other models to see a variation in the lensing magnification effect due to 
possible changes in the shape of the redshift distribution. 

The upper panels of Figure |l| show the skewness parameter 53 (#) evaluated without taking the magnification effects into 
account. As the figure clearly shows, the skewness is very sensitive to both f2 m and the mean redshift of sources, but insensitive 
to the shape of the redshift distribution (comparing model A with B). 

The lower panels of Figure |l| show the skewness correction due to the lensing magnification effects. We took a(z) — 
which gives the strongest estimate of the lensing magnification effect as will be discussed in §4. The nonlinear result deviates 
from the linear one on a sub-degree scale and the nonlinearity reduces the amplitude of the magnification effect as it has 
stronger influence on V K (9) than on («g) mag . Between the scales displayed, the skewness correction is almost constant, the 
variation is less than 1.4. The lensing magnification effect becomes stronger as the mean redshift becomes lower and as the 
distribution becomes broader. This redshift distribution dependence is, at lest qualitatively, similar to that on the source 
clustering effect, which may be due to the similarity in their phenomena (Hamana et al. 2000). The most important point 
found in Figure nl is, however, that the correction term is small, 3% correction at largest. 




4 DISCUSSION AND CONCLUSION 

We have examined the lensing magnification effects on the convergence skewness using the nonlinear semi-analytic approach. 
Numerical computations were done only for the case of a(z) — 0. Does this choice have a special meaning ? So far, little is 
known about the number counts of distant galaxies with measured redshifts, in particular at z > 0.2 where most of the source 
galaxies are located. We can, however, put a constraint on the range of the "effective value" of a (under the assumption of 
no strong redshift evolution in a) as follows: (i) a > by definition, (ii) The observational facts that the number density of 
the distant galaxies behind the lensing galaxy clusters in optical or near IR bands is smaller than that measured in the field 
region (e.g., Fort et al. 1997). This indicates that a < 1 (comes from eq. (Q) with fi > 1 within cluster region). Therefore, it 
may be said that the possible range of the "effective value" of a is < a < 1. Within this range, a — gives the upper limit 
of the lensing magnification effect on the convergence skewness. We may, therefore, conclude that the lensing magnification 
has no significant effect on the convergence skewness, its amplitude could be 3% at most. 

In the above discussion, we implicitly assumed that the source galaxies are selected by a flux threshold. However, one 
can adopt the selection with a surface brightness criterion. In this case, the flux magnification has no influence on the cosmic 
shear measurements, but the area magnification still exists. Therefore, the lensing magnification effect in this case can be 
estimated by setting a = 0. Thus, the effect is stronger for the surface brightness selection than for the flux selection. 



© .... RAS, MNRAS 000, §-0 



6 Takashi Hamana 




Figure 1. Upper panels: Skewness parameter S3 of the lensing convergence evaluated without taking the lensing magnification effect 
into account. The dot-dashed line is the quasi-linear perturbation theory computation for A model. The other curves represent the 
semi-analytic computations with the nonlinear matter bispectrum fitting formula (vWHSCBOl). Lower panels: Skewness correction term 
due to the lensing magnification effects (for a case of a(z) = 0). The dot-dashed line is the linear theory computation for A model. The 
other lines represent the semi-analytic nonlinear computations. 



Finally, we notice a limitation of our calculation. We used the weak lensing approximation to derive equation (|3|). This 
approximation works well except for very rare directions such like the core of clusters of galaxies where the lensing surface 
mass density is very high. Such regions where a strong lensing event can take place are very small, 9 ^ 0.5 arcmin, at largest. 
Therefore, it is expected that the strong lensing may change the above results on scales below 1 arcmin. 
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